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A multiple-graphene-layer (MGL) structure with a stack of GLs and a highly conducting bottom 
GL on SiC substrate pumped by optical radiation is considered as an active region of terahertz 
(THz) and far infrared (FIR) lasers with external metal mirrors. The dynamic conductivity of the 
MGL structure is calculated as a function of the signal frequency, the number of GLs, and the 
optical pumping intensity. The utilization of optically pumped MGL structures might provide the 
achievement of lasing with the frequencies of about 1 THz at room temperature due to a high 
efficiency. 



I. INTRODUCTION 



Since the first experimental demonstrations of the non- 
trivial properties of graphene (see, for instance, Refs. [l{ 
2] and numerous subsequent publications), due to the 
gapless energy spectrum of graphene, the latter can be 
used as an active region in the terahertz (THz) and far 
infrared (FIR) lasers with optical or injection pump- 
ing 0, 0] • The optical pumping with the photon energy 
Ml leads to the generation of electrons and holes with 
the energy £o = Ml/ 2. Since the interaction of elec- 
trons and holes with optical phonons is characterized by 
fairly short time To, the photogenerated electrons and 
holes quickly emit cascades of N optical phonons, where 
N = [eo/hdJo] and [X] means the integer part of X. Due 
to relatively high energy of optical phonons in graphene 
(hu)Q ~ 200 meV), the number TV can vary from zero 
(pumping by CO2 lasers) to a few units (pumping by 
quantum cascade lasers or semiconductor injection diode 
lasers). Thus, the photogenerated electrons and holes 
populate the low energy regions of the conduction and 
valence bands of graphene [|[ (see also Refs. H, @]). As 
a result, the Fermi levels of electrons and holes are sepa- 
rated and shifted from the Dirac point to the conduction 
and valence bands, respectively. This corresponds to the 
population inversion for the interband transitions with 
absorption or emission of photons with relatively low en- 
ergies: hw < 2sf, where ef is the quasi- Fermi energy 
of the electron and hole distributions, and leads to the 
negative contribution of the interband transitions to the 
real part of the dynamic conductivity Re er w (which in- 
cludes the contributions of both interband and intraband 
transitions). If Re a u < 0, the stimulated emission of 
photons with the relatively low energy hw (in the THz 
or FIR range) is possible. The intraband contribution is 



primarily associated with the Drude mechanism of the 
photon absorption. As shown previously the realization 
of the condition of lasing is feasible [1, @] ■ Some possi- 
ble schemes of the graphene lasers (with metal waveguide 
structure or external metal mirrors) utilizing the above 
mechanism of optical pumping were considered recently 
(see, for instance, Ref. @)- 

In particular, the laser structure evaluated in Ref. 0], 
includes two graphene layers. Relatively weak absorption 
of optical radiation with the efficiency f3 = ire 2 /he ~ 
0.023, where e is the electron charge, H is the reduced 
Planck constant, and c is the speed of light, per one pas- 
sage necessitates the use of rather strong optical pump- 
ing. This drawback can be eliminated in the struc- 
tures with multiple graphene layers (MGLs). Such epi- 
taxial MGL structures including up to 100 very per- 
fect graphene layers with a high electron mobility (/1 ~ 
250,000 cm 2 /sV) preserving up to the room tempera- 
tures and, hence, with a long momentum relaxation time 
were recently fabricated using the thermal decomposi- 
tion from a 4H-SiC substrate [1, The incident op- 
tical radiation can be almost totally absorbed in these 
MGL structures providing enhanced pumping efficiency. 
A long momentum relaxation time (up to 20 ps |10j ) im- 
plies that the intraband (Drude) absorption in the THz 
and FIR ranges can be weak. The optically pumped 
MGL structures in question with the momentum relax- 
ation time about several ps at the room temperatures 
might be ideal active media for interband THz lasers. 
The situation, however, is complicated by the presence 
of highly conducting bottom GL (and, hence, absorbing 
THz and FIR radiation due to the intraband processes) 
near the the interface with SiC as a result of charge trans- 
fer from SiC. 

In this paper, we analyze the operation of THz lasers 
utilizing optically pumped MGL structures by calculat- 
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FIG. 1: (a) Schematic view of a laser with a MGL struc- 
ture, (b) Occupied (by electrons) and vacant states in differ- 
ent GLs under optical pumping. Arrows show transitions re- 
lated to interband emission and intraband absorption of THz 
photons with energy ftio (interband transitions related to op- 
tical pumping as well the processes of intraband relaxation of 
the photogenerated electrons and holes are not shown). 



cupied. Taking into account that elevated electron and 
hole densities (at elevated temperatures and sufficiently 
strong optical pumping considered in the following), the 
electron and hole distributions in the range of energies 
e <C fiujo in the k—th graphene layer (1 < k < K) can 
be described by the Fermi functions with the quasi-Fermi 

energies 6p [see Fig. 1(b)]. The case of relatively weak 
pair collisions of the photogenerated electrons and holes 
(due to their low densities) in which the energy distri- 
butions deviate from the Fermi distributions was studied 
in Ref. [|[. Using the Falkovsky-Varlamov formula [ll[ 
for the dynamic conductivity of an MGL structure gen- 
eralized for nonequilibrium electron- hole systems Q , one 
can obtain 
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ing their characteristics as function of the number of 
graphene layers K and optical pumping intensity Iq and 
demonstrate the feasibility of realization of such lasers 
operating at the frequencies from about 1 THz to several 
THz at room temperatures. 



II. MODEL 

We consider a laser structure with an MGL struc- 
ture on a SiC substrate serving as its active region. Al- 
though the active media under consideration can be sup- 
plemented by different resonant cavities, for dcfimtcncss 
we address to the MGL structure placed between the 
highly reflecting metal mirrors (made of Al, Au, or Ag) as 
shown in Fig. 1(a). It is assumed that the MGL structure 
under consideration comprises K upper GLs (to which 
we refer to just as GLs) and a highly conducting bot- 
tom layer with a Fermi energy of electrons £p, which is 
rather large: Ep ~ 400 meV @. Apart from this, we 
briefly compare the MGL structure in question with that 
in which the bottom GL is absent. The latter structure 
can be fabricated by an additional peeling of upper GLs 
and depositing them on a Si substrate. 

Due to the phogeneration of electrons and holes with 
the energy £o = M7/2 followed by their cooling associ- 
ated with the cascade emission of optical phonons, low 
energy states near the bottom of the conduction band 
and the top of the valence band can be essentially oc- 



for the GLs with 1 < k < K. Here t b and r are the elec- 
tron and hole momentum relaxation times in the bottom 
and other GLs, respectively, T is the electron and hole 
temperature, and k B is the Boltzmann constant. The 
first and second terms in the right-hand sides of Eqs. (1) 
and (2) correspond to the interband and intraband tran- 
sitions, respectively, which are schematically shown by 
arrows in Fig. 1(b). For simplicity we shall disregard the 
variation of the electron and hole densities in the bottom 
GL under the optical pumping, so that the electron-hole 
system in this GL is assumed to be close to equilibrium 
with the Fermi energy Ep determined by the interaction 
with the SiC substrate. 

The quasi-Fermi energies in the GLs with k > 1 are 
mainly determined by the electron (hole) density in this 
layer i.e, oc and, therefore, by the rate of 

Ik) 

photogeneration G n by the optical radiation (incident 
and reflected from the mirror) at the k—th. GL plane. 
Using Eq. (2) for fuo = Ml, we obtain 
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Here 1^ is the intensity (power density) of the optical 

pumping radiation at the k— th GL. At Ml > 2ep (for 

all GLs), Eq. (3) yields ~ 01^ /Ml. Considering 
the attenuation of the optical pumping radiation due to 
its absorption in each GL, one can obtain 



P) K - k +(i-fi B Y0--m 



(4) 



Here Iq is the intensity of incident pumping radiation and 
Pb = (47t/c)<Tq. The latter quantity accounts for the 
absorption of optical pumping radiation in the bottom 
layer. 

A relationship between and Gq ^ is determined by 

(k) 

the recombination mechanisms. We assume that e F oc 



[Gq' 1 ] 7 , where 7 is a phenomenological parameter, 
this case, 
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where e f = e F ^ is the quasi-Fermi energy in the top- 
most GL. 



III. MGL STRUCTURE NET DYNAMIC 
CONDUCTIVITY 

Taking into account that the thickness of the MGL 
structure is small in comparison with the wavelength of 
THz/FIR radiation, the generation and absorption of the 
latter is determined by the real part of the net dynamic 
conductivity: 



K 



Re cr w = Re a„ + Re <r, 



(fc) 



(6) 
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Taking into account that in the frequency range under 
consideration hoj -C £ F) so that one can neglect the first 
term in the right-hand side of Eq. (1) responsible for 
the interband absorption in the bottom GL, and using 
Eqs. (1), (2), and (6), we arrive at 
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FIG. 2: Frequency dependences of the real part of dynamic 
conductivity Re a w normalized by quantity e 2 /Ah for MGL 
structures with different number of GLs K at modest pump- 
ing (ej = 15 meV.) 



the interband transitions. When the latter is negative, 
i.e., when the interband emission prevails the interband 
absorption, the quantity Re as a function of ui ex- 
hibits a minimum. At a strong optical pumping when 

(k) 

the quantities e F are sufficiently large, Re a u < in 
this minimum as well as in a certain range of frequencies 



< U! < UJ n 



Here uj r 



are the frequen- 



cies at which Re = 0; they are determined by r B , 
t, and Ep' (i.e., by the intensity of the incident optical 
pumping radiation). 

Since in the MGL structure in question Ep ^> k B T, 
considering such frequencies that uj 2 t b , uj 2 t 2 ^> 1, one 
can reduce Eq. (7) to the following: 
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Under sufficiently strong optical pumping when e f 
Hui, k B T and, consequently, tanh[(fio> 



-1, setting J2k=i £ 
Eq. (8) one obtains 



- (fc) ~ K*e%, where K 
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The first two terms in the right-hand side of Eq. (7) de- 
scribe the intraband (Drude) absorption of THz radia- 
tion in all GLs, whereas the third term is associated with 



IV. FREQUENCY CHARACTERISTICS OF 
DYNAMIC CONDUCTIVITY 



Figures 2 and 3 show the frequency dependences of 
Re <7 W normalized by e 2 /Ah calculated for MGL struc- 
tures with different K at different values of £p (i.e., differ- 
ent optical pumping intensities) using Eqs. (1), (2), (5), 
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FIG. 3: Frequency dependences of the real part of dynamic 
conductivity Re normalized by quantity e 2 /Ah for MGL 
structures with different number of GLs K at £§• = 30 and 
50 meV. The inset on upper panel shows how the GL popu- 
lation varies with the GL index k, whereas the inset on lower 
panel demonstrates the dependences but in a wider range of 
frequencies. 
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FIG. 5: Real part of the normalized dynamic conductivity as a 
function of the number of GLs for different frequencies. The 
line with markers corresponds to a MGL structure without 
the bottom GL. 



positive at w > uj max , where Lo max is rather large: more 
than 10 THz. The MGL structures with larger number 
K of GLs at stronger pumping exhibit smaller uj m m and 
deeper minima Re a^. This is confirmed also by Fig. 5. 



V. ROLE OF THE BOTTOM GL 
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FIG. 4: Dependences of u>min on number of GLs K at differ- 
ent e F . 



and (6) or Eqs. (5) and (7). We set ef = 400 meV Q, 
fcl = 920 meV, T = 300 K, tb = 1 ps, r = 10 ps, and 
7 = 1/4. As seen from Figs. 2 and 3, Re a u can be nega- 
tive in the frequency range u > uj m i n with ui m i n decreas- 
ing with increasing quasi-Fermi energy in the top- 
most GL, i.e., with increasing optical pumping intensity 
(see below). In the MGL structures with K — 50 — 100 
at ep — 30 — 50 meV, one has cj min /2Tr > 1 THz (see 
Figs. 3 and 4) . As follows from the inset on upper panel 

(k) 

in Fig. 3, the quantities e F are not too small in compar- 
ison with e 1 ^, even in GLs with the indices k -C K, i.e., 
in GLs near the MGL structure bottom. This implies 
that the pumping of such near bottom GLs is effective 
even in the MGL structures with K ~ 100. The quantity 
Re cr w as a function of oj exhibits a minimum (see the in- 
set on lower panel in Fig. 3). The sign of Re cr w becomes 



The highly conducting bottom GL, whose existence is 
associated with the intrinsic features of the MGL struc- 
ture growth, plays the negative role. This is mainly be- 
cause it results in a marked absorption (due to the Drude 
mechanism) of the THz radiation emitted by other GLs. 
Since such an absorption increases with decreasing fre- 
quency, the achievement of the negative dynamic con- 
ductivity in the MGL structures with shorter momentum 
relaxation times tb and r even at higher frequencies is 
significantly complicated by the Drude absorption. This 
is demonstrated by Fig. 6 (upper panel). In contrast, a 
decrease in the Fermi energy ej? and the electron density 
in the bottom GL, of course, might significantly promote 
the achievement of negative dynamic conductivity in a 
wide frequency range, particularly, at relatively low fre- 
quencies (compare the frequency dependences on upper 
and lower panels in Fig. 6). In this regard, the MGL 
structures with GLs exhibiting a long relaxation time r 
(like that found in Ref. [T(J) but with a lowered electron 
density in the bottom GL or without the bottom layer ap- 
pears to be much more preferable. Such MGL structures 
can be fabricated using chemical/mechanical reactions 
and transferred substrate techniques (chemically etching 
the substrate and the highly conducting bottom GL [12| 
or mechanically peeling the upper GLs, then transferring 
the upper portion of the MGL structure on a Si or equiva- 
lent transparent substrate) . The calculation of Re a u for 
this MGL structure can be carried out by omitting the 
term Re a u in Eq. (6). The pertinent results are shown 
in Fig. 5 (see the marked line) and Fig. 7. Here as in 
Fig. 3, we assumed that that hw = 920 meV, T = 300 K, 
tb = 1 ps, and r = 10 ps. 

The obtained frequency dependences are qualitatively 
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FIG. 6: Real part of the normalized dynamic conductivity 
versus frequency calculated for MGL structures with tb = 
0.1 ps and different r and ep. 
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FIG. 7: Real part of the normalized dynamic conductivity 
versus frequency calculated for MGL structures with different 
numbers of GLs and without bottom GL. 



similar to those shown in Figs. 2 and 3. However the 
dependences for the MGL structures without the bottom 
layer exhibit a marked shift toward lower frequencies. In 
particular, as follows from Figs. 5 and 7, Re can be 
negative even at ui/2ir < 1 THz (at choosen values of 

el)- 



VI. CONDITION OF LASING 

To achieve lasing in the MGL structures under consid- 
eration, the following condition should be satisfied Q: 



8tt 

— Re a u 

c 



E 2 > (l- ri )E 2 + (l-r 2 )E 2 - 



(a/R) 2 E 2 +E 2 . 

(10) 

Here, E t , E\, and E% are maximum amplitudes of the 
THz electric field E = E(z) at the MGL structure (placed 
at the distance t from the the bottom mirror, where t is 
the thickness of the substrate), and near the pertinent 
mirror, respectively, Eg = (asns/2) J Q E 2 dz, as and 
ns are the absorption coefficient of THz radiation in the 
substrate (SiC or Si) and real part of its refraction index, 
r\ and T2 are the reflection coefficients of THz radiation 
from the mirrors, and a/R is the ration of the diame- 
ters of the output hole a and the mirror R. In deriving 
inequality (10), we neglected the finiteness of the MGL 
thickness (in comparison with t and the THz wavelength) 
and disregarded the diffraction losses. For simplicity, one 
can set E^Ef ~ E 2 and E 2 S - (ta s n s /2)E%, disregard- 
ing, in particular, partial reflection of THz radiation from 
the bottom GL. In this case, inequality (10) can be pre- 
sented as 



8tt 



Re a u | > (1 - n) + (1 - r 2 ) + (a/R) 2 



ta s n s /2 = L. 

(11) 

Assuming that r% = r-i = 0.99, (a/R) = 0.1, as — 
2 — 4 cm -1 , ns — 3 [13], and t = 50 fim, for L one 
obtains L = 0.06 — 0.09. However, as follows from Fig. 5, 
the quantity (87r/c)|Re cr w | for a MGL structure with 
K = 100 at uj/2it = 1.5 THz at e T F = 30 meV is about 
12/3 ~ 0.275. At oj/2tt = 1.0 THz but for the struc- 
ture without the bottom GL (see the line with markers 
in Fig. 5), one obtains (87r/c)|Re er w | ~ 0.345. These val- 
ues of (87r/c)|Re cr^l well exceed the above value of L (in 
contrast with the structures with two GLs Q, for which 
the minimization of the losses is crucial) . At the elevated 
frequencies, the ratio (87r/c)|Re cr M \/L can be even much 
larger. 

Considering that the electron (hole) density E T in the 
topmost GL 



2 

^h 2 



dpp 



1 + exp[(v F p - e^)/k B T] 



at e T F = 30 meV at T = 300 K, we obtain E T ~ 
2 x 10 11 cm~ 2 . Such a value (and higher) of the photo- 
generated electron and hole density is achievable exper- 
imentally (see, for instance, Ref. [1J])- Assuming that 
K = 100, Hfl = 920 meV, and the recombination time 
t r ~ 20 ps at T = 300 K and E T ~ 2 x 10 11 cm" 2 [H[, 
for the pertinent optical pumping power we obtain 1^ ~ 
6.4 x 10 4 W/cm 2 . One needs to point out that this value 
of 1^ is much larger than the threshold of lasing at a 
certain frequency (uj m in < u < w max ). At T = 100 K, 
the recombination time (due to optical phonon emission) 
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is much longer [15[, hence, the electron and hole densi- 
ties in question can be achieved at much weaker optical 
pumping. 

In the regime sufficiently beyond the threshold of lasing 
when the stimulated radiative recombination becomes 
dominant, the pumping efficiency is determined just by 
the ratio of the energy of the emitted THz photons 
ftto and the energy of optical photons Ml: r\ = uj/fl. 
However, the nonradiative recombination mechanisms 
can markedly decrease rj. In the regime in question, 
the maximum output THz power can be estimated as 
maxf w ~ 7r R 2 (ui/tt)lQ. For example, at ft£l — 920 meV, 
fkj/2-K = 5.9 meV(w/27r ~ 1.5 THz), 2R = 0.1 cm, and 
In = 3 x 10 4 W/cm 2 , one obtains max P u ~ 1.5 W. 



VII. CONCLUSIONS 

We studied real part of the dynamic conductivity of 
a multiple-graphene-layer (MGL) structure with a stack 
of GLs and a highly conducting bottom GL on SiC sub- 
strate pumped by optical radiation. It was shown that 
the negative dynamic conductivity in the MGL structures 
under consideration with sufficiently large number K of 
perfect upper GLs can be achieved even at room temper- 
ature provided the optical pumping is sufficiently strong. 
Due to large K, the absolute value of Re<7 w in its mini- 
mum can significantly exceed the characteristic value of 
conductivity e 2 /4H. Thus, the MGL structures can serve 
active media of THz lasers. This can markedly liberate 
the requirement for the quality of the THz laser resonant 
cavity. An increase in tb and r promotes widening of the 
frequency range where Rec^ < 0, particularly, at the low 
end of this range. This opens up the prospects of THz 
lasing with lu/2tt ~ 1 THz even at room temperature. 
The main obstacle appears to be the necessity of suffi- 
ciently long relaxation times tb and r in GLs with rather 
high electron and hole densities: E > 10 11 cm -2 . Since 
the electron-hole recombination at the temperatures and 
densities under consideration might be attributed to the 
optical phonon emission [1 61 ] . the optical pumping inten- 
sity required for THz lasing can be markedly lowered (by 
orders of magnitudes) with decreasing temperature. The 
conditions of lasing in the MGL structures without the 
bottom GL are particularly liberal. 
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Appendix A. Recombination 

The rate of radiative recombination in the degenerate 
electron-hole system in the topmost GL due to spontan- 
ious emission of photons at e% 3> ksT can be calculated 
using the following formula [g, EE[ 



11 ' frr / = S&r <* (4) 3 . (Al) 



2v r 



Trh 3 



3TTh 3 V 3 F 



Here v r = yfaz(&e 2 /3hc){yF/c) 2 VF and ae is the dielectric 
constant pf — s'f/vf- 

The rate of the electron-hole recombination associated 
with emission of o ptic al phonons can be described the 
following equation 14j : 



Rph cx / dpp(hu! /v F -p) 



1+exp 



v F p ■ 
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-exp 
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oc exp 
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k R T 
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Equalizing R^ and the rate of generation of electrons 
and holes by the optical pumping radiation Gq and the 
recombination rate, and considering Eqs. (Al) and (A2), 
one can find that £^ oc (G^) 1 / 3 and oc ln(G^), respec- 
tively. The radiative interband transitions stimulated by 
the thermalphotons can also contribute to the recombi- 
nation rate [a, [l5[ as well as the processes of electron- hole 
interaction in the presence of disorder. These processes 
provide different dependences of the recombination rate 
on the quasi- Fermi energy. Due to this, calculating the 
dependence of the latter in different GLs on the optical 
pumping intensity, we put, for definiteness, oc (G^) 7 
with 7 = 1/4. Such a dependence is less steep than that 
for the spontaneous radiative recombination and some- 
what steeper than in the case of optical phonon recom- 
bination. The variation of parameter 7, leads to some 

(k) 

change in the e y F — k dependence, but it should not af- 
fect the main obtained results. 
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